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Introduction 


There is a close relation between the concepts of entropy and viscosity, 
associated with systems of conservation laws. It is well known, for example, 
that vanishing viscosity weak solutions for such systems must satisfy the 
entropy inequality across their discontinuity, and that the converse holds, at 
least in the small - in the large for scalar problems; both are being used to 
identity the so-called "physically relevant" solution of such systems, e.g., 
[ 7 ]. 

In this paper we amplify a certain aspect of this relation, with regard 
to conservative difference schemes 

( 1 . 1 ) 

v^(t+k) = v^(t) - ^[h(vy_p^j(t),***,Vy^(t)) - h(vy_p(t),*‘*,v^_j(t))] 
serving as consistent approximations to the scalar conservation law 
(1.2) (x,t) +-||- (u(x,t)) = 0. 

To make our point, consider a difference scheme which is known to satisfy the 
entropy inequality; roughly speaking, this should indicate according to the 
above, the existence of certain amount of numerical viscosity present in such 
a scheme. It is therefore plausible to assert that other schemes, containing 
more numerical viscosity, will also have to satisfy the entropy inequality. 
After putting these terms in a more precise framework, we show the validity of 
the above assertion subject to the technical assumptions listed below. Thus, 
we prove the entropy inequality by means of comparison. 
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In [12], Osher introduced, for the method of lines, a class of E 
schemes which were shown to converge to the physically relevant solution. 
Making use of the terminology just introduced, the so-called E schemes can 
be identified as exactly those having no less numerical viscosity than that of 
Godunov. Since the latter is known to satisfy the entropy inequality, we are 
able to extend Osher ^s ideas to the fully discrete case, as a special case of 
the above assertion. This is carried out in Section 5, paving the way for the 
proof of the more general assertion in Section 6. Prior to that, we give in 
Sections 3 and 4, a brief discussion on the entropy inequality in relation to 
the all Important Godunov and Lax-Frledrichs schemes. 


2 • Preliminaries 

We consider difference schemes 


(2.1a) Vy(t-Hc) = H(vy_p(t),-«*,v^p(t); f,X) 


which admit a conservative form 


(2.1b) 

H(Vy_p.-*-.v^p; f.A) = v^ - ^[h(v^_p^j,*--,v^) - h(v^p,-‘-,v^_j)]. 

and are serving as consistent approximations to the scalar conservation law 
(2.2) (x,t) + (u(x,t)) = 0. 


Here, v^(t) H v(x^,t) is denoting the approximation value at the gridpoint 



(x^ = vAx,t), k and Ax are respectively, the temporal and spatial meshslze 
such that the mesh ratio X = k/Ax is being kept fixed, and p, a natural 
number. Finally, h^ly,^= h(v^_^ 1 1 * * * * ^ Lipschitz continuous 

numerical flux consistent with the differential one, h(w,w,»**,w) = f(w); for 
the sake of simplifying the notations, its possible dependence on f and X 
is being suppressed. 

We begin by putting the scheme (2.1) in an Increment form; using the 
difference operator - w^, we set for v^^j^ * 


(2.3a) 


(2.3b) 


V 2 




= X 


^v+1 ~ ^v+ 


adding and subtracting f^ to the RHS of (2.1b) and making use of (2.3), 
(2.1a) reads 


(2.4) Vy(t+k) = Vy(t) + c;;;^ I /2 I /2 (t) - (t). 


Next, we denote 


(2.5) 


= C,.. 1 , + C.., 1 , = X 


V I/2 ~ v+ I/2 ^ V2 


f V fy+l ~ 1/y 


Av 


V+I /2 


Af If 

Noting the identity - C^i. = X ■ ^ '2 

y V+I /2 V1-1/2 ^^v+1/2 


, the Incremental coefficients 


'vf I/2 


equal 


C* = T (Q 1^-+ ^ )• 

2 ^^v+ 1/2 


(2.6) 
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Inserted into (2.4), our scheme then recast Into the form 
(2.7) 

Vy(t-Hc) = v^(t) [f(v^+l(t)) - f(vy_i(t))] +Y 

which reveals the role Q plays as the numerical viscosity coefficient. We 
will therefore use Q as a measurement of the amount of viscosity present in 
such a scheme. 

REMARK ; In the case of 3-polnt schemes, p “ 1, this measure of viscosity is 
rather general in the sense that such schemes are completely determined by 
their coefficient of numerical viscosity, e.g., [10]. We do not claim the 
generality for (2p+l) -point schemes p > 1: this definition of numerical 

viscosity is in fact 3-polnt oriented, as we shall see in a more precise form, 
later on. 

Let TV[v(t)] = |v^j^(t) - v^(t)| denote the total variation of the 

computed solution at time t; we then have the following 


T-Arnma 2.1 

The scheme (2.1) has a total variation non-increasing provided its 
numerical viscosity coefficient, Q^^l/ > satisfies 


( 2 . 8 ) 


Af 

A^ 


y+V2 

v+V? 


< Q 


v+ lA 


< 1 . 


Proof ; The inequalities (2.8), expressed in terms of the Incremental 
coefficients in (2.6), are translated into 
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(2.9) 






> 0 , 


1 - C 


V+I/ 2 " °v+l/2 


> 0 . 


k straightforward calculation, based on the Incremental form (2.4) and the 
Inequalities (2.9), shows the non-increase In total variation, 

TV[v(t+k)] < TV[v(t)], see [5]. 

Lemma 2.1 implies, in particular, the convergence of the scheme (2.1), 
provided its numerical viscosity coefficient meets the requirement (2.8): one 

can select a boundedly a.e. converging subsequence, v (t;Ax^), such that its 


limit v(x,t) 


lim V (t;Ax^) satisfies (2.2) in the weak sense. 


x=vAx' ,Ax'-»^0 

[9]^^^. Weak solution of (2.2) however, are not unique. The 
lower-bound on the LHS of (2.8) requiring that much of viscosity for 
convergence to ^ limit weak solution, does not guarantee this weak solution to 
be the physically relevant one: it is well-known for example, that the 3- 

I V? 1 

point Courant-Isaacson-Rees scheme where ^ j , may admit 

limit weak solutions violating the physically relevant entropy condition, 
e-g*, [5], [12]; thus, a greater amount of viscosity is required for the 
entropy condition to hold. In the next section we discuss Godunov's scheme 
which turns out to play a central role in determining that additional required 
amount • 

We note in passing, the fundamentally different role played by the upper 
bound on the numerical viscosity, appearing on the RHS of (2.8). It is 
related to the hyperbolic nature of the approximated equation (2.2), as it 
amounts to the CFL-like condition, see (2.5), 


( 1 ) 


We consider compactly supported initial data; a further L ^bound. 


j 7 ^ j . r i ■■■ 7 , CO ^ 

derived below, is required for the more general Initial data in L 0 L fl BV. 
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(2.10) ^|(f(Vy) - h^l/ 2 ) + " Vl/ 2^1 ^ ” ""vl* 

which usually results In placing a limitation on the mesh ratio, A, being used 
(recall that h(***) may depend on A as well). A stricter CFL condition of 
this type was Introduced in [9]. In particular, the numerical flux of a 
difference scheme satisfying (2.10) admits the consistency relation 

(2.11) *^(%-p+l»***»’^v-l*''»”»’^v+2*'**»^v4p) “ 

Such essentially 3-point schemes include, beside the standard 3~polnt schemes, 
several of the recently constructed second-order accurate converging schemes, 
e.g., [5], 18]. 

Finally, we would like to point out that by halving the CFL number, one 
obtains a maximum principle; that is 


Lennaa 2.2 



satisfying 


( 2 . 12 ) 




Then, the following maximum principle 


(2.13) 


lnf[v (t)] < V (t-Hc) < sup[v ,(t)] 
P P 


holds 
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Proof ; The Incremental coefficients In (2.6) do not exceed a value of 


0 < 1, - (q 1/ + ^ ^ ) < •^•2Q 1# < 4 

v+1/2 2 '‘^vf V2 ^ ^ ^V2 2 


Making use of the Incremental form of the scheme, see (2.4), 


,(t+k) - C^i/2V^^j(t) + (l - 


and noting the convexity of the combination on the RHS, (2.13) follows. 


3* The Entropy Condition and Godunov's Scheme 

The building block In Godunov's scheme, [4] , Is the solution of the 
Rlemann problem. Let u^(x/t; denote the similarity solution 

of the Rlemann problem (2.2) subject to Initial conditioning 

u(x,t=0) = y + 1 + sgn(x) ^ Godunov's scheme Is deter- 

L lert 4 rignt 

mined by 


(3.1a) 

where 

(3.1b) 

(3.1c) 


v^(t+k) 


H^(v. 


v-r 


v+1 


) - 2 ’ 


^v-1/2 = 


1 ^~x/2 « 

^ / u (x/k;vy_j,vy)dx 

1 R 

/ u (x/k;Vy,Vy^j)dx 


Assume the CFL condition 
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(3>2a) A»max |a(u)| < ^ * a(u) = f(u) 

u 


holds* The RHS 
see Figure 3-1, 

of (3.1b-c) can be evaluated from the 
giving 

integral form of (2.2), 

(3.2b) 

- 2X (f^,- 


(3.2c) 



where 



(3.3) 


stands for the 
(3.2b-c), (3.1a) 

numerical flux of Godunovas scheme; 
takes the desired conservative form 

Indeed, by averaging 

(3.4) 

V^(t+k) = V^ - ^lly^) 

• 


^v- 1/2 V2 1/2 



( ^ Ax/2 ■^“1 ^ Ax/2 < 1 ^ Ax/2 -| ?• Ax/2 <■ 

Figure 3-1 
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Conslder a pair of scalar functions (u(w),F(w)) such that 


(3.5a) 


U(w)f(w) = F(w), 


U(w) > 0; 


the entropy condition for a physically relevant solution of (2.2), u = u(x,t), 
requires the following entropy Inequality 


(3.5b) 


U(u) +|^F(u) < 0 (weakly) 


to hold for all entropy pairs related through (3.5a). Recalling (3.1b-c), 
Jensen's Inequality and the Integral form of (3.5b) yield 


Ax/ 2 


(3.6a) U(vy^l, ] < / U(u*^(x/k;v^_j,v^))dx < U(v^) - 2 X(f(v^ - Fv-Vo^ 


(3.6b) U(v^_^ly^) < ^/^^U(u^(x/k;vy,v^j))dx < U(v^) + 2X(F(Vy) - ) 


where 


(3.7) ^S+1/2 ^ = F(u^(0''’;vy,v^j)) 


is Godunov's numerical entropy flux, consistent with the differential one, 
F^(w,w) = F(w). Averaging (3.6a-b) we find on account of (3«la), that 

Godunov's scheme is consistent with the differential entropy inequality 
(3.5b), 


(3.8) U(vy(t-Hc)) < 


U(v‘ 


,.U 


) + U(v 
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We now summarize what we have shown In the following 


Lemma 3*1« 

Assume the CFL condition 


(3.9) 

X’max |a(u) 1 < -^ 


u 

holds. For 


(3.10a) 

“ G 

%-l/2“ % - - ^v-l) “ Qv-1/2^''v-1/2 

(3.10b) 

Vv+l/2= - ^(^v+1 - fv) + Q\h-1/2^^v+1/2 


we have the following entropy Inequalities 

(3.11a) - 2 A(f(v^) - ) 

(3.11b) ^ + 2 A(f(v^) - 

Proof : Inserting the definition of the numerical viscosity coefficient in 

(2*5) , one reads (3«10a-b) from (3m2b-c)# The conclusion appears in (3«6a-b)m 

REMARK; We have shown that Godunov's scheme satisfies the entropy inequality 

(3 ,8) by averaging (3.11a-b), while assuming the CFL condition (3.2a), 

A*max|a(u)| < ; the latter was required in order to guarantee that waves 
u 

issued from the two opposite faces of the v^-cell, do not interact. In the 
scalar case, an entropy solution is known to exist whether or not these waves 
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interact# Hence, (3.8) follows from the integral form of (3.5b) applied over 

the whole v^-cell (rather than — as we have done — over its left and right 

halves), provided the relaxed CFL condition X»max|a(u)| < 1 holds, thus 

u 

preventing these waves from reaching the cell's other faces. The reason for 
our introduction of v^(t4k) as the average of ^^ 1/2 » each of 

which satisfies the entropy inequality (3.11a-b), will prove itself essential, 
however, in studying E schemes in Section 5 below. We note that the so 

Introduced averaging is nothing else but a restatement of the following 
identity whose verification is left to the reader 


(3.12) = 




In closing this section, we would like to point out the following 
geometric interpretation of the numerical viscosity coefficient associated 

Q 

with Godunov's scheme, * integrating (2.2) over the left half of 

v^^l-cell, see Figure 3-1, we find 


v^+V = - 2X[f 


v+1 “ "v+ 


vJ 


while Integration over the right half of v^-cell yields, as before. 


V^+l/ = + 2X(f^ - ]. 


Subtracting the second from the first, we have 


V^+I/ - %+l - % - + ^v+1 " 2^5+1/ ) = (1 - 2Q^1, )'(vy+i- v^), 
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Thus, 1 - 2Q^_j_Xyr^ gives us the compression ratio 

(*WV2‘ 


4. Lax-Frledrlchs Scheme and Its Entropy Satisfying Hodlflcatlon 
Lax-Frledrlchs scheme [2], [6], given by 


(4.1) 


- H (v„.pV^,v^jiX) 5 -j 


V+1‘ 


’v-1- 


has the most allowable numerical viscosity under the total variation non- 

LF 

Increasing requirement (2.8), Q lyr^ = 1. A. Harten has observed [private 
communication] that the scheme coincides with that of Godunov, when the latter 
Is applied over a staggered grid, see Figure 4-1, 


Ax 


v^(t-Hc)=-^ / u^(x/k;Vy_j,Vy^j)dx = H^^(v^_j,v^,v^j;X), 

-Ax 


provided the CFL condition 


(4.2) 


X*max|a(u)| < 1 
u 


Is met 
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Vy(t+k) 



Integrating the differential entropy inequality (3.5b) over the same domain, 
we end up with its discrete version 


U(v^(t-Hc)) 


U(Vy+l(t)) + U(v^_j(t)) 


- 2 ~ 


after little rearrangement, it can he put into the more standard form, compare 
(3.8), 


(4.3a) 


< U(v„) - 


where 


(4.3b) f!Ji/2= = 2 


, ) + F(v^) , 


is LF numerical entropy flux, cons;lstent with the differential one 
F^^(w,w) = F(w). 
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We note that LF scheme does not admit a simple averaging of the type 
Introduced above for Godunov's scheme. Instead, one might consider the 
following modification 


(4.4) 


V . . + 2v.. + V. 


v^(t-Hc) - ^ f 


The so-modified scheme has half the numerical viscosity of LF- 
M 1 

scheme, \j^ = V2 » and can be rewritten in the desired averaged form 


(4.5a) v/t*) . 


where 


(4.5b) 


+ ^v-1 




I/2 = = 2^ " “ Vl) 


(4.5c) 


+ - +1 


'v+l + 




The new scheme introduced, (4 #4), can also be viewed as a two-cell averaging 
of two non-interacting Riemann problems, see Figure 3-1, 


^ . Ax/2 p 

Vy(t-Hc) = h“(v^_j,v^,v^j) [ / u (x/k;v^_j,v^)dx 

-Ax/ 2 


j. 

+ \ u (x/k;Vy,v^j)dx] 


provided the CFL limitation 
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(4.6a) X*max|a(u)| < ^ 

u 

is met. Integrating the entropy inequality (3.5b) over the same two-cell 
domain, the scheme is found to satisfy that entropy inequality is its standard 
discrete version 


(4.6b) 


U(v^(t-Hc)) < U(v^(t)) - - Fv-1/2^ 


with a numerical entropy flux 


(4.6c) F^i/ = 


+ F(v„) 


- - (■X'-v+i) -“<%» 


M 

consistent with the differential one, F (w,w) = F(w). 

In analogy with Lemma 3, we are now ready to state 


Lemma 4*1 

Assume the CFL condition 


(4.7) 


X*raax|a(u)| < -i 


u 


holds « 


I2L ^v±l/2» -S4yen_:^^ 


%_!, = Vy - - fy_i) - Qy_l/ Avy_i/ 


(4.8b) 


%+lA = - ^(^v+1 - ^v) + » 


we have the following entropy inequalities 
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(4.9a) 

D(v;.l/^) < 0(»„) - 2A(F(»„) - 


(4.9b) 



Proof; 

Recalling that 0^1/ * y » Ihe RHS of (4,8a) and (4.8b) 

coincide 

with 

“*^(''v'''\n-l’''v+l’^') reapectlvelyi 

applying 


for the latter the LF entropy inequality as quoted in (2.3a-b), one reads the 
conclusion (4.9a-b). 


REMARK ; In [7], P. D. Lax gave a direct proof of the entropy inequality 
(4.3), for LF scheme approximating arbitrary system of conservation laws. (In 
comparison, the arguments used in the above scalar analysis, requires the 
existence of an entropy satisfying Rlemann solution in the larse .) Since the 
modified scheme is nothing else but an average of two LF-solvers, Lax's result 
goes over in this case; that is, for arbitrary system of conservation laws, 
both LF and the modified scheme, satisfy the entropy Inequality for all 
entropy pairs associated with the differential system^^\ As much as we are 
aware, these are the only two known examples satisfying the entropy condition 
in such generality. 


5. The Entropy Condi ton and E Schemes 

In this section we study difference schemes containing no less numerical 

Q 

viscosity than that of Godunov, Q > Q • Such E schemes — after Osher [12] 


^^^For the exact CFL limitation in this case, see [7]« 
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— are shown to converge to the unique physically relevant solution of (2.2), 
provided the CFL limitation 


(5.1) 




) + 



Is met. Ideally, one would like to allow the relaxed CFL limitation (2.10) to 
be used; the reason for introducing the stricter (5.1) — half the usual CFL 
number — stems from the fact that we were unable to rewrite LF scheme in the 
desirable averaged form as discussed In Section 4. We note that (5.1) takes 
the equivalent form 


(5.1') 



< 


1 

T * 


which, in the case of Godunov's scheme, amounts to preventing waves 
interaction. As before, such a CFL limitation yields, in particular, the 
consistency relation (2.11), characterizing essentially 3-point schemes. 


Theorem 5.1 

An E scheme converges to the physically relevant solution of (2.2), 
under the CFL restriction (5.1). 


Proof ; Convergence was established in Lemma 2.1, since an E scheme is 
necessarily total-variation non-increasing 


X 


Af 

Av 


v+ V? 
V+I/2 


< Q°+V 2 ' V1/2 " 


2 


< 1 . 


We turn to examine the entropy inequality. We attach the superscript G, M, 
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and E to distinguish between Godunov's scheme (3.1), the modified scheme 
(4.4) and the E scheme under consideration (2.7) 

Vy(t+k) = Vy(t) -j (f(v^j(t)) - f(vy_j(t))) +-| (t))). 

We rewrite the latter in the averaged form 


(5.2) 



where 



(5.3a) 

E- 

%- 1/2 

“ - X(f^ - f^_i) - Q^l/^Av^l/^ 

(5.3b) 

E+ 

^ v + 1/2 

= Vy - X(f^^j - f Av^ 


Recall the corresponding averaging forms for Godunov's scheme, see (3.10), 
(5.4a) = v^ - X(f^ - f^_i) - Q®_ 1 / 2 Av^_ 1 /^ 

(5.4b) = Vy - X(f^^ - f^) + 

and that for the modified scheme, see (4.8), 

(5.5a) v“: - X(£„ - f^.j) - ^ iv„. 

'w 1/2 ■ ’'v ■ ■ *v) + i l"'v+l/2 ■ 


(5.5b) 
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Accordlng to our assumption 


(5-6) (1 - ) I . 


0< ®V±14‘ 1- 


Multiply (5.4a) by 6y_ , (5.5a) by (1 - 0^ ) and add to find that 

(5.3a) amounts to 


VvZl/2 “ ®v-1/2^vIi/ 2+ - ®v.ly^)^vIl/2 5 


similarly, multiplying (5.4b) by , (5.5b) by (1 - 0^^ ) and adding, 

we end up with (5.3b) having the form 


E+ n G+ ^ „ 0 \ M+ 

Vv+1/2 = V 1/2 1/2+ (1 - V1/2 ^''v+1/2* 


Averaging the last two equalities, (5.2) reads 


(5.7) 


V (t+k) 
V 


V-V 2 G- ~ ^v- 1 / 2 ^ M- 

2 V 2 ^ 2 ^v- V 2 


9 , (1 - 0 1 ) 

v+ 1/2 G+ v+ V2 M+ 

+ V , + V , 

2 v+ 1/2 2 v+ 


Thus, we see that every E scheme can be written as a convex combination of 
one-sided averaged Rlemann solutions * 

Let (u(w),F(w)) be an entropy pair associated with (2.2). By the 
convexity of U, (5.7) Implies 


(5.8) D(v„(ttk)) < -!^0(v°: “(''vI \ ) 
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Next, we invoke the entropy Inequalities concluded in Leimnas 3.1 and 4,1 


U(v^l/J < U(v^) ” 2 X(f(v^) - F®_i/ ) 


uCvS^l/,) < U(vy) + 2X(F(uy) - ) 


U(v“:i/^) < U(v^) - 2 X(f(v^) - f“_ 1 /^) 
U(v^l/, ) < U(v^) + 2 X(f(v^) - 


When Inserted into (5.8), we end up with the desired entropy inequality 


(5.9a) 


uK(t-Hc)) < U(v^) - >^(F^ 1 / - F^_i/ ) 


with a numerical entropy flux 


(5.9b) 


^v+ V2 ~ 1/2 ^''+1/2'^ ^ ^ “ ®V+ I/2 ) ^ v+ 1/2 


consistent with the differential one, F (• • • ,w,w, • • • ) = F(w). 


REMARKS; 

(i) An explicit formula for Godunov's numerical flux, 

“ Min[sgn(Vy^j - v^)f(v)], < v < was given in [12]; 

here = Min/Max(Vy,v^j ). Hence, an equivalent characterization 

for E schemes, requiring 


sgn(Vy^l - Vy)[h^l, - f(v)] < 0, 


min . . max 

Vv+1/2 < V < 
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shows that a 3-polnt monotone scheme Is an E scheme. Unfortunately, E 

schemes like monotone ones, are at most first-order accurate [12]. 

(11) We have seen that E schemes satisfy the entropy Inequality (5.9) 

for all entropy functions, U(*); their corresponding ntmerlcal entropy fluxes 

are given as convex combinations of two numerical fluxes associated with 

monotone schemes -- Godunov and the modified LF scheme (4.4). Hence, an Lj- 

\u 

convergence rate estimate of order (4x) ^ follows along the lines of [9, 

Theorem IV] 

lv(*,t) - u(*,t)l , < Hv(*,t=0) - u(*,t“0)II +K*(tAx) ^2. 

L L 

Considerations of the constant coefficients case shows this L^-estlmate to be 
sharp, e.g., [11, Sections 9 and 10]. 

(Ill) As an Immediate corollary from 'Hieorem 5.1 we obtain verification of 
the following "folklore" result. 


Corollary 5.2 

A conservative difference scheme with a non-vanlshlng nimerlcal 

viscosity, 0 < Q . < Q,,, 1# (^) ^ X » Is converging to the unique entropy 

mill v+ •*-/2 2 

solution for sufficiently small mesh ratio, X. 

Such non-vanishing viscosity schemes were specifically "tailored”, for 
example, in [5, Section 5]* Here we note, that the CFL-like restriction on 
the mesh ratio, X, depends heavily on the behavior of the flux, f, near the 
sonic points • 
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6. Ruaerlcal Viscosity and the Entropy Condition 

In this section we would like to systemlze the kind of arguments 
Introduced above, emphasizing those essential ingredients which prevail in the 
more general context. 

We consider a general conservative scheme which we rewrite in the 
averaged form, compare (2.7), 


( 6 . 1 ) 

v^(t+k) 


[v (t)-X(f -f )-Q , Av , ]+[v (t)-X(f -f )-K) ^1 Av 1 ] 

V V v-1 V- V2 V2 V ^+1 V V2 ^ V2 

» 


the entropy condition follows by constructing a consistent discrete entropy 
inequality for each of the averaged terras on the RHS (6.1), thus opening the 
door for showing the former by means of comparison . For that purpose, pick a 
3-polnt entropy condition satisfying scheme 


(6.2) v (t+k) » h(v ,v ,v ;f,x) = v - x(h “ h i, ) > 

v -1 V v +1 V '-VH- V2 V 2 


such that the following holds: 


Assumption : The numerical flux is independent of the mesh ratio X. 


The plausibility of the above assumption stems from the fact that the Rlemann 


problem admits similarity solution u^(x,t;u 




) 


left right 
R 

u hence all difference approximations based on 

Riemann solvers must satisfy such a requirement; this is not the case, for 
example, with LF scheme (4.1), where we were forced to consider instead its 


modification (4.4). 
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The reason for Introducing the last assumption Is becoming clear upon 
writing 


(6.3) 


- X(f - f ) - Q Av . ] 

V V -1 ~V- 1/2 V- I/2 


[v - X(f - f ) + Q 1 , Av 1 , ] 
^ V v+1 V ~v+ 1/2 ^ V2 ^ 


where, see (2.5), 


(6.4) 


Q 1, = ^ 

~v+ 1/2 


f + f - 2h 
V v+l 


v+ V 2 


• Av 


vf I/2 


depends linearly on 
abbreviated as before 
as 


X; hence, the two averaged terras on the RHS of 
by v^^ and — can be equivalently 






(6.3) — 

expressed 


each of which satisfies the entropy Inequality, provided the CFL limitation Is 
being halved # Termwise comparison of the averaged forms, (6.1) and (6#3), 
shows their difference only In the numerical viscosity coefficients; assuming 
^v+l /2 vary between two coefficients of numerical viscosity associated 

with entropy satisfying schemes, we are able to represent (6*1) as a convex 
combination of the latter. The discrete entropy Inequality follows for the 
corresponding convex combination of entropy fluxes. 


We have shown 
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Theorem 6.1 

Consider the difference scheme ( 6 » 1 ) and assume that the CFL condition 


IQ 


V+I/2 




^ - ^v +1 ~ 


Av, 


\H-l/2 


<1 


holds . Then, the scheme satisfies the entropy condition, provided we can 
find another entropy satisfying difference approximation with less numerical 
dissipation , , 

V2 ^ V2 * 


The corresponding numerical entropy flux Is given by 


I/2- Q 


_ r ^+V2 l., . r ^‘/2 


Q 1 - Q 1 

v+ 1/2 ~\H- 1/2 -I M 


^^^One may assume. Instead, iQ^+l/^l 4 lQ^+l/2l, ^+1/2 denoting the 

numerical viscosity coefficient of a difference scheme admitting the desired 
entropy satisfying averaged form. 
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